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Abstract 



O The prescription of Silva to derive superpotential equations from variational deriva- 

J tives rather than from Lagrangian densities is apphed to theories of gravity derived from 

i_^ Lovelock Lagrangians in the Palatini representation. Spacetimes are without torsion and 

T— H isolated sources of gravity are minimally coupled. On a closed boundary of spacetime, the 

> 

On metric is given and the connection coefficients are those of ChristofFel. We derive equations 

^^ for the superpotentials in these conditions. The equations are easily integrated and we 
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^ of gravity. Finally, we find a simple expression for the superpotential of Einstein- Gauss- 

es 

Bonnet theories with an anti-de Sitter background: it is minus the KBL superpotential, 
confirming, as it should, the calculation of the total mass-energy of spacetime at spatial 
infinity by Deser and Tekin. 



give the general expression for all superpotentials associated with Lovelock Lagrangians. 
We find, in particular, that in Einstein's theory, in any number of dimensions, the super- 
potential, valid at spatial and at null infinity, is that of Katz, Bicak and Lynden-Bell, the 
KBL superpotential. We also give explicitly the superpotential for Gauss-Bonnet theories 
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1 Introduction 

(i) A view on superpotentials 

In electromagnetism, one of Maxwell's equations relates the stationary electric field E to 
the density of charges pe'. V ■ E = ^irpe- From this follows that the total charge Q responsible 
for the field is equal to the flux of E through a closed surface surrounding the sources. E 
is the electric force acting on a unit test charge q. In Newton's theory of gravitation the 
gravitational field G is related to the density of matter pm in a similar way: V • G = —Airpm- 
Thus the flux of — G through a closed surface surrounding the source responsible for G is 
equal to the total mass M. G is the gravitational force acting on a unit test mass m. In 
Einstein's theory of gravitation things get slightly more complicated. First there is a change 
of meaning: total mass is now total mass-energy Mc^. Second, given a localized source of 
gravity, part of the total mass-energy is in the gravitational field itself though its density 
is not defined nor is it possible in general^ to even disentangle the total gravitational field 
energy from the total mass-energy of spacetime. Nevertheless, an isolated amount of matter 
together with its gravitational field appears from a great distance as a point-like source of 
gravitation, possibly spinning, and like in Newton's theory and in electromagnetism its mass- 
energy is also equal to a flux across a closed surface at spatial infinity. However a surface 
element in spacetime is a 2-index anti-symmetric tensor so that instead of a vector (E or G) 
we have a two index antisymmetric tensor whose flux across the surface at infinity equals 
mass-energy. The tensor is commonly called the superpotential. It is worth noting that the 
flux includes the energy of the sources [20]. 

Unfortunately the superpotential is not as well defined as in classical field theory. Given a 
Lagrangian from which Einstein's equations are derived, there exists a unique superpotential 
associated with that Lagrangian^. The trouble is that the Lagrangian density itself is only 
defined up to a divergence and is thus not unique. That is the reason why efforts have 
been made to obtain the superpotential not from the Lagrangian but from the variational 
derivatives of the Lagrangian^. These are insensitive to additional divergences. 



^See however [21] and [22]. 

■^This is not shown here. It can be derived from the Belinfante-Rosenfeld identities, see [20], without 
invoking boundary conditions or field equations. 

^And not from field equations as is sometimes stated, see for instance [13]. 



(a) Superpotentials derived from variational derivatives 

Silva [30] suggested a method to obtain superpotentials from variational derivatives. Pre- 
vious works, in particular those of Joe Rosen [29] and of Anderson and Torre [2] aimed at 
obtaining conserved quantities from the field equations. Here the emphasis is on superpo- 
tentials and there is no need for field equations. The fluxes we are after are in addition not 
necessarily conserved like, for instance, the Bondi mass [3] at null infinity. In the last decades 
relativists have become interested in spacetimes with more than 4 dimensions or which far 
from the sources are not flat; they become, in particular, anti-de Sitter^. The proper way to 
deal with this situation as well as in keeping the ability to calculate, for instance mass-energy, 
in asymptotically flat spacetimes in spherical coordinates, is by introducing a background 
metric^ [17]. Ferraris, Francaviglia and Raiteri [13], who generalized, to some extent, Silva's 
idea were concerned with the variations of conserved quantities. 

(Hi) What is done in this paper 

Silva's Lagrangian approach is summarized in Section 2. Julia and Silva [18] applied the 
method to Einstein's relativistic theory of gravitation and found that given the asymptotic 
metric components at spatial or at null infinity, the unique answer is the KBL superpotential 
[20]. Silva's prescription applies to variational derivatives that depend at most on first order 
derivatives of the field components and whose Lie derivatives with respect to some arbitrary 
vector field depend also at most on first order derivatives of this vector field. These conditions 
are met in a GL(D, R) formulation of Einstein's theory. Ferraris, Francaviglia and Raiteri 
[13] extended Silva's prescription to Lagrangians that contain higher order derivatives of the 
fields but, in the Lie derivatives, do not contain higher orders of derivatives of the vector 
field. Their work renders Silva's prescription applicable to Einstein's ordinary equations in 
terms of the metric components and the superpotential we derived that way is again the 
KBL one. In Section 3 we extend Silva's prescription to the case in which Lie derivatives 
contain second order derivations of the vector field. This makes the method applicable to 
variational derivatives of Lovelock Lagrangians® [24] (Section 4) in the Palatini representation 



''The mass-energy in such spacetimes have been calculated in a series of fine works by Deser and collabo- 
rators [1], [8], [9] to which we shall come back later. See also the paper by Petrov [27]. 

^Elsewhere backgrounds appears in the guise of "counter-terms" like in [15] or a "regularization" procedure 
as in [23] . A detailed discussion on the background formalism is given in [20] . See also an interesting discussion 
in [18]. 

®The method is equally applicable to more general higher order derivative Lagrangians. 



in which the metric and the symmetric connection coefficients are the independent fields. We 
find the equation of the superpotential and its solution. The superpotential of Einstein's 
gravity theory in N dimensions is, as expected, the KBL superpotential (Section 5). The 
superpotential for Gauss-Bonnet variational derivatives, derived in Section 6, is far more 
complicated that the KBL superpotential. The superpotential of Einstein- Gauss-Bonnet 
theories is a linear combination of both. In particular, the superpotential with an anti-de 
Sitter background is proportional to the KBL superpotential. The one associated with the 
spherically symmetric solution given in [1] is minus the KBL superpotential and on a sphere 
at spatial infinity this coincides with the superpotential found by Deser and Tekin [9]. 

(iv) A word about identities 

It may be noticed that nowhere are field equations used. From the beginning to the end we 
play with identities. That is not altogether surprising. Globally conserved or non-conserved 
quantities are constants or functions that appear in the boundary conditions. 

It is also interesting to note that the derivation of the superpotential from a Lagrangian 
via Noether identities, like it is done in [19] or [20], is very different from the derivation from 
variational derivatives. We have not yet tried to relate the two methods; the relation appears, 
to us at least, somewhat mysterious. 



(v) What has this to do with conservation laws? 



This paper is far removed from conservation law considerations and is more of a mathe- 
matical character. It is therefore perhaps useful to reconnect superpotentials to conservation 
laws before we move on. The superpotential is an antisymmetric tensor density, linear in 
some arbitrary displacement vector say'^ ^ , defined on the boundary of the domain, usually 
at spatial or null infinity where spacetime identifies with the background that does not have 
to be flat. The relation to "conservation laws" is as follows, at least in Einstein's theory of 
gravitation. The ordinary divergence of that antisymmetric tensor density is a divergence- 
less vector density, say J^. If one uses Einstein's field equations, and this is one of the rare 
occasions where they are used in this paper, it was found® [20] that J^ = \/—gTaS,^ plus 
terms which in the linear approximation are negligible. If, in particular, the background has 



^For a definition of notations see Section 2 below. 
*When spacetime is flat far from the sources. 



a Killing field of time translations like in a flat spacetime, then one clearly sees that the 
total flux of the conserved current is related to mass-energy conservation. However, if the 
background possesses no Killing field then JP is a conserved vector with no obvious physical 
meaning. So much for conservation laws and the relevance of superpotentials. 

2 The Silva prescription 

(i) Elements 



In what follows we are given a set of tensors on a curved spacetime with components 



"^ and a Lagrangian density^ 



The particularity of C is that variational derivatives, 

d£ „ dC 



(2.11 



C 



Da- 



(2.2) 



dvA d{D„yA) 

contain only first order derivatives of yA- Notice that partial derivatives with respect to yA 
are at DuyA constant and vice versa. The variation of C with respect to yA can thus be 
written in this form: 



6£ = C^SyA + Do 



dC 



d{D,yA) 



^VA 



C^5yA + d^ 



dC 



d{D„yA 



-^VA 



(2.3) 



The spacetime has here the passive role of a background so that 5D = D6. The Lie derivative 
of a Lagrangian density^*^, 

£^C = D^iCn = dA^n- (2.4) 

Thus, if (5 = £^, (2.3) minus (2.4) must be equal to zero: 

dC 



£^C - 5,(£r) = C''£^yA + d. 



d{D„yA) 



£m - ^^" 



0. 



(2.5) 



^In Einstein's theory of gravity which is in 4 dimensions, Greek indices go from to 3. Latin indices from 
1 to 3. The signature of the metric p^^ is —2 and g is its determinant. Covariant derivatives are indicated 
by a D; partial derivatives by a d. The permutation symbol in 4 dimensions is e^,^p„ with eoi23 ~ 1 and in 3 



dimensions by ekim with ei23 — 1. The 4- volume element dV^ 



31 c^i/po- 



da;" A dx'' A dx"^ , a surface element in 



4 dimensions dS^,^ — ^e^vpadx'' A dx" and in 3 dimensions dSt = \^okimdx^ A dx"^ = ^ekimdx' A da;'". The 
Sections on Einstein's gravity and Gauss-Bonnet gravity is in A^-dimensions to which these definitions extend 
in an obvious way. In those Sections, Greek indices run from to A'' — 1 and the signature of the metric is 
-(iV-2). 

^°One should pay attention to the fact that we deal with vector densities and that ordinary divergences of 
such vectors as well as of anti-symmetric tensor densities are scalar or vector densities. 



(2.5) is commonly referred to as Noether's identity. Next we assume that £^yA has this form 

£iyA = ^Axt + KxDp^^- (2.6) 

This is the case for tensor fields. The A's depend on i/a and D\yA and are tensors. Inserting 

(2.6) into C £^yA, we see from (2.5) that this expression, which we designate by X, is 
necessarily of the following form: 

X ^ C^'l^yA = C^Aax^^ + ^""KxDpe = d. [j^^ + U^'Dp^) ■ (2-7) 

The J's and C/'s are components of tensor densities. A derivation by part of the second term 
after the first equality sign leads to what has been called since the late 1940's generalized 
Bianchi identities because this is what they are if yA is the gravitational field: 

C^Aax - dp{C^A%) = 0. (2.8) 

All this has been derived in innumerable papers but has been re-derived here to clarify 
our notations. 

(ii) The "cascade identities" of Julia and Silva 
In their paper of 1998 [17], £/^ is replaced by 

e = Co^e, (2.9) 

where ■^0(2;^) is regarded as "fixed" for a moment and e{x^) is an arbitrary scalar function. 
Inserting (2.9) into (2.7) and expanding in terms of e, 9pe and Dp„e, one obtains for the 
left-hand side of (2.7), an expression of this form 

X = Xoe + W(ldpe where W^ = C^A%Co- (2.10) 

Xq is X in which ^ has been replaced by ^q. In the right-hand side of (2.7) we have a 
divergence of 

J^C^ + U^'D,C^ = JSe + U^'dpe with J^ ^ Jx^o+U^Dp^o and U^' ^ t/^eo- (2.11) 

But following (2.7), (2.10) is equal to the divergence of (2.11) and since ^q is arbitrary the 
identity also holds if we replace ^jy by any ^ (we remove then the indice 0). Thus, 

Xe + W^dpe = da {re + U^^dpe) = d^J^e + (J" + D^U'^p) d^e + U^P^^Dp^e. (2.12) 



e being arbitrary, one can identify the coefficients of e, dpe and Dp„e of both sides: 

X = dpJi" , WP = JP + d^U^P and U^P''^ = 0. (2.13) 

The interesting point is that U^'^ is anti-symmetrical, U^" = f/''"^', and the identities may be 
rewritten 

X = dpJP , JP = WP + daUP''. (2.14) 

If C = 0, then X = 0, W^ = and the vector density JP is the divergence of an 
antisymmetric tensor density: 

JP = d^UP'' , (only if £^ = 0). (2.15) 

(2.15) will not be used later. 

(in) Functional equation for a superpotential 

Now first for the motivation. Consider the integral, over a volume V with boundary S, 
of JP as given by (2.14): 

/ jPdVp = I WPdVp + / UP^dSpa. (2.16) 

Jv Jv Js 

In particular, let the volume be a spacelike hypersurface, say x^ = 0. (2.16) can then be 

written 

0j3 f wO j3_ , / TTOk 



rd'x = / Wd'^x + * U'^'dSk- (2.17) 

V Jv Js 

Let further ^ be a timelike vector which on the boundary, at spatial or null infinity, is 

associated with time translations in a fiat background. The timelike field exists in the whole 

spacetime (no black holes). In coordinates in which ^ = {1,0,0,0} the volume integral 

becomes 

/ J^d^x = [ C^A%d^x + / U'^'^dSk. (2.18) 

Jv Jv Js 

In general relativity, (2.18) is the Hamiltonian: 

/ J^d^x = / 7=5(^11 - KT^)d^x + / U'^'^dSk with GP = RP - I 5PR, (2.19) 

Jv Jv Js 

Rpa is the Ricci tensor and R the scalar curvature while T^ is the energy momentum of 

the source of gravity and k is the usual coupling constant^^. When the field equations are 



^ About K see (4.7) below and the corresponding footnote. 



satisfied, the Hamiltonian is equal to the total mass-energy: 



,2 _ J, TTpcr^a — X TjOkjc^ i 1,, ;r rA 



Mc^ = <i UP^'dSp^ = <i U'^^'dSk , (only if £^ = 0). (2.20) 

Js Js 

JJP'^ is thus the superpotential by definition. 

The Hamiltonian field equations are obtained by applying the variational principle to 
the Hamiltonian. Regge and Teitelboim [28] brought attention to the fact that if boundary 
conditions are given to begin with, the variation of the Hamiltonian should have no boundary 
terms like in classical mechanics. Hamiltonians have no time derivatives. Silva's prescription 
is a covariant expression of that remark. So here is how it goes. 

Consider the variation of (2.16) due to arbitrary variations of yA', it can be written 

dWP ^ 

a(zw)'^^- 

f2.21) 



/ Jf'dVp] = I ^^6yAdVp + J) iy/P" + bU'"') dSp^ where W^" 



The condition that no boundary term appear in (2.21) provides a functional differential 
equation for the superpotential: 

bUP" = -W^'' so that bil Jf'dVA = / ^^^byAdVp. (2.22) 

Notice that the surface integral in (2.21) and equation (2.22) are only correct if W^^'^' = 0. 
That this is indeed true can be seen as follows. Consider the variational derivatives of X, 
which according to (2.13) on the left, must be equal to zero: 



6X dX 



Sy, ay, °'(to)=" '""" ^ = ^'^^^e + C-A%D,(\ (2.23) 

Now in (2.23) replace S,'^ by ^Qe, expand as a polynomial in e, dpe and Dpo-e and equate 
the factors of e, dpe and Dpa^e to zero. Then remove the indice because ^g is arbitrary. The 
factor of Dp^e is W^^'^'. Thus W"' is indeed antisymmetrical. 

Next suppose that W'' is a solution of equation (2.22) on the left. Then, U^'' + CP'^, 
where C^^^^x"^) is an arbitrary function independent of yA, is also a solution of (2.22). When 
the superpotential relates, in particular to mass-energy, C"^ defines its "zero" point. As it 



was done in [20] and also in [18], we take C"^ = —U^'^ which is the superpotential of the 
background. The superpotential of the background is obtained from U^'^ by equating to zero 
the source terms-'^^. The final solution is then given by U^^ = V^ — V^" . If, however, we use 





^^For instance, in the Schwarzschild spacetime we take m = in whichever coordinates the metric is written. 
Notice that if spacetimes are asymptotically flat, V'^dSpa = while if they are anti-de Sitter, V" dSpa = oo 

but in both cases V" dSpc, is bounded as can be seen in the example considered at the end of the last Section. 



a background from the very beginning, as we have done here and shall do in our examples, 



then U^^ = and U^^ = JJP" . So much about the works of Silva and that of Julia and 


Silva. For completeness let us mention a work of Fatibene, Ferraris and Francaviglia [12] 
which deals with relative conservation laws. 

3 Higher order derivatives of ^ 

The method described so far is not applicable to the Palatini formulation of General Relativ- 
ity because the Lie derivatives of the connection coefficients contain second order derivatives 
of ^ . So, we first generalize the prescription of Silva to the case where Lie derivatives contain 
second order derivatives of ^ . 

(i) Basic elements 

Equations (2.1) to (2.5) remain valid here but instead of (2.6) we assume that 

£^yA = ^Axe + ^''AxDpe + ^'I^D^^p„^e => Kx = Ky (3-1) 

(2.7) has now this form 

X = C^£^yA = d^ [r,e + Ul'D^e + VT^"^(M.)^') • (3.2) 

Set again, like in (2.9), ^ = ^jye. The left hand side of (3.2) can then be written like this, 
thanks to (3.1), 

X = Xoe + W^dpe + yP'^Dp^e with < ^ C^Kx^o + '^Kx^aCo and Y^''' ^ ^^A^^Co • 

(3.3) 
As before, the indice indicates that ^o is there instead of ^^. Thus (3.2) takes the following 
form 



Xoe + Wlldpe + yP'^Dp^e = 8^ ^J^e + t/o^^a^e + V^^'^'D^^ej where V^^'"' = V^'^" . (3.4) 

Since (3.4) holds for any ^q, it also holds for any ^^ and 

Xe + W^dpe + V^Dp^e = {d„r)e + {JP + D„U''P)dp€ 

+ (i/(p-) + DxV'^nDp^e + V^P'^Dy^p^e. (3.5) 

(a) Cascade equations 



In local Minkowski coordinates, Dxpa^ = d\pcj^ and its factor, V^^P''' in the right-hand 
side of (3.5), must be zero because there is no similar factor on the left-hand side and e is 
arbitrary. Thus, we must have both 

y\pcT ^ y\ap ^^^ y{\pa) ^ 1 ^y\pa _^ ypa\ _^ yaXp^^ ^ q. (3.6) 

Multiplying (3.6) by Dxp„e, we get 

(^Xpa ^ ypaX ^ yaXp^^ ^^^^^ ^ ^ (37^ 

With (3.6), one finds that (3.7) can equally be written as follows: 

V^'^Dxp^e = -lV''''''R\p,dpe, (3.8) 

and with (3.8), we can now rewrite (3.5) like this 

X€ + WPdpe + YP''Dp^e={d„r)e + {JP + D„U''P-lV^"'^R\^^)dp€ + {U'^fi''HDxV^P'')Dp^€. 

(3.9) 
This identity holds for any e. Thus the factors of e, dp€ and Dpue from both sides of the 
equality must be identical, i.e. 

X = d^r , WP = JP + D^tj^'P - IVP^^R^^^ and YP-" = U^p^'^ + DxV^P''. (3.10) 

These cascade equations are similar to (2.13). We shall now reduce this set to one that looks 
exactly like (2.13). 

(in) An equivalent set of identities 

We start from WP as given in (3.10) in which we replace U'^p by {—U^'^^ + U'^P"'^): 

WP = JP- d^tj^P"^ + D^U^'^P'^ - lVP^^R\^^. (3.11) 

Then we replace D„U^^P' by its value deduced from YP^ , given in the right-hand side of 
(3.10), and we get for (3.11) 

WP = JP - d^U^P''^ + D„Y"P - [d^xV^p'' + IVP^^RP^^^ . (3.12) 

However, thanks to the symmetries of V P"^ we readily find that 

D^xV^P^ + lVP^^R\^^ = -|a. (DaV^'^'I^) . (3.13) 

10 



Thus instead of (3.12) we may also write 

WP = JP - d^UP" + D^YP"" where UP" = U^P''^ - IDxV^P"^^ = -U"P. (3.14) 

The tensor [/'"^ is anti-symmetrical. Thus (3.14) provides identities similar to (2.13): 

X = dpJP , JP = *WP + d^UP'' where *WP = WP - D^YP" . (3.15) 

and the corresponding generalized Bianchi identities are 

C^Aax - DpiC^A'^^) + Dp^iC^'AZ) = 0. (3.16) 

(iv) More identities 

Consider now (3.3) which we rewrite like this: 
X = Xo + W^dpe + YP'^Dp^e = Xq + Wppe + YP^dp^e where WP = WP - YP^p^,. (3.17) 
Since X is a divergence we have 

'■^ "-^ dJJ^)=0. (3.18) 



SyA dyA \d{D„yA) 

Inserting (3.17) into (3.18), identifying the coefficients of e,dpe,dpo-€ and dxpo-e to zero, re- 
moving the in dice 0, because the resulting identities are valid for any ^'^, one obtains four 
identities, two of which do not involve X: 

f)VP^ flWP f)Yf"^ 

-^5yA = W^"''^ => WP" ^ T^ ,SyA and Y^^P"^ = => Y^P" ^ —- -5yA. 

dyA d{DayA) o{DxyA) 

(3.19) 

We notice that Y P" has the same symmetries as V P" . Therefore (3.13) holds for Y P" 

except that in (3.17) we have ordinary derivatives as if we were in a flat space and therefore 

daxiY^n = -§5. {dxY^P"^^) . (3.20) 

(v) The superpotential equation 

We now consider the variation of JP as given in (3.15), using WP as defined in (3.17) 

SJP = 6{WP - D^YP") + d„5UP" = 5WP - dJYP" + dJUP". (3.21) 
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Notice that W^ like Y^^^ contains at most first order derivatives of yA- Therefore, we may 
write (3.21), making use of (3.17) as follows: 



5WP 
OVA 



dYP" ( dyp" \ 



(3.22) 



With (3.19), the sum of the first two terms in parenthesis may be written 



and therefore bJ^ can also be written 

bJP = - — by A + d^m^P"^ - dxij^''") + bVP^V (3.24) 

OVA 
But using (3.20) and the definition of W in (3.17) we find that 

d^W^P'^^ - d^x^^'"') = d^^^P"^ + iBxiJ^P''^^)^ (3.25) 

WP" has been defined in (2.21). So 

^jp = ^^^byA + d^W^P''"' + |L>a(^''"'1^) + ^UP% (3.26) 

byA 

But, see (3.15) and (3.17), 

bWP b*WP 
WP = WP-YP''TP^^ = *WP + d^YP^ so that ^ = ^V— • (3-27) 

P- byA byA 

Thus the differential equation for the superpotential is now 

bUP'' = -W^P"^ - iDxiY^"^^^) and b( f JPdVp^ = f ^^byAdVp. (3.28) 

Putting together the various elements defined along the way, we have an equation for the 
superpotential - the first equality in (3.28) - that looks like this: 

f)W^ ( rlY^^P \ 

d{D^]yA) \d{D^]yA) J 



in which 



W = £^A^;^^^ + 2£^A^'"^Z)^e^ and y^'^ = r^A^'^^C^. (3.30) 



Equations (3.29) and (3.30), together with appropriate boundary conditions, are all we need 
to calculate the superpotential in the following example. 
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4 Application to variational derivatives of Lovelock Lagrangians 
in the Palatini representation 

(i) Lagrangian, variational derivatives and boundary conditions 

In the Palatini formulation we take the inverse metric components g^^'^ and the connection 
coefficients T'*' = T\ as independent field components. The curvature tensor does not 
have all the symmetries of the Riemannian one. To avoid any confusion we shall denote 
the curvature tensor, like in Eisenhart [10], by B\ and reserve the usual notation of the 
curvature tensor R\p^ for when the connection coefficients are Christoffel symbols. B\ 
is antisymmetrical in the last two indices only, B\ = —B\^. Additional symmetry 
properties, similar to those of the Riemann curvature tensor, are: 

B\upa) = , B\pa:r) = « and D^B^, = DpB,„ - D„B,p where B„, = B\^, / B,„. 

(4.1) 
Here a semi-column means covariant differentiation. The curvature tensor itself is given by 

B\p, = 2{d^pT\^^ + T\pT\^^). (4.2) 

In accordance with [20], we introduce a background whose metric components are 5 . Instead 
of r's we introduce 

A\ ^ T\, - f\, = \ g'^^'iDpg^, + D^g^p - D^gp,). (4.3) 

The bars like on D^s refer to the background. A is a tensor. The curvature tensor can be 
written in terms of these A's: 

S V = mP^\]. + ^\[p^%) + ^V- (4-4) 

As boundary conditions we want to impose the value of the metric components on the 
boundary S, i.e. 

9,iu\s=gpu- (4-5) 

We shall also demand that on the boundary the T's be Christoffel symbols i.e. that 

Dxgns = ^ and thus B\p„\s = R\p^. (4.6) 
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The action of the gravitational field is of the form 

^ = ^ / l^{%v. B\,)<f^ where k = '^ll^^L. (4.7) 

£ is a Lovelock Lagrangian [24] , N is the dimension of spacetime. The coupling constant k 
is normalized like in [1], except for a factor 2: S'iv-2 is the surface of a sphere of dimension 
N — 2 and Gn the gravitational coupling constant. ^^ 

Contrary to our assumption in Sections 2 and 3 that spacetime has a passive role here 
it is part of the game: the metric components and the connection coefficients are the field 
components. Nonetheless, Section 3 is applicable to variational derivatives of Lovelock La- 
grangians because: (1) the connection coefficients enter into tensorial combinations, (2) the 
variational derivatives as well as Lie derivatives of the connection coefficients are tensors as 
well and (3) on the boundary the connection coefficients are Christoffel symbols. 

The variation of the action has this form 

5 A = JiCf.uSg'"' + £/"(5A^^Jd^x + f V^'dS^. (4.8) 

C^u and C^'^ are variational derivatives. In applying the variational principle to the action 
of the gravitational field plus its sources, {5 A + 5^sourccs = 0), £^1^ is defined by the sources 
but for minimally coupled matter, C^'^ = 0. This last equation is linear and homogeneous in 
Dxg^'^ and contains no derivatives of B\ . This is the particularity of Lovelock Lagrangians. 
As a consequence, and as shown by Exirifart and Sheikh- Jabbari [11], Dxg^'^ = will always 
be a solution of this equation. We shall, of course, not use that solution except to notice that 
it is in accordance with the bounday condition (4.6). In addition to C^'^ = 0, the variational 
principle would impose some boundary condition i.e., for isolated sources, 

VdSf, = 0. (4.9) 

Adding a divergence to the Lagrangian may change this condition but, contrary to what 
happens in the Hamiltonian formalism, there is no way to get rid of a condition like (4.9) by 
adding a divergence to C. The condition (4.9) may be satisfied if (4.5), (4.6) hold. Otherwise 
(4.9) must at least be compatible with (4.5), (4.6). 

(ii) The superpotential 



^For Af = 4, & = 47r and k. = ^^ in which G is Newton's gravitational constant. 
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The Lie derivatives £^yA are as follows, the second equality is a repeat in terms of A's 
like in (3.1) where indices A are surrounded by curly brackets: 



The only non-zero A's appearing in W^ and y"^, see (3.30), are 



{lJ.v}ri 



Al^^r.^_2gV(Ps';^} and Af;l^^l = USiiS:; + 6^^6i^)5l 



Introducing these A's into W^ and Y^'^ , we obtain 



(4.10) 
(4.11) 



(4.12) 



WP = 2{-C{^^ + C.P''D„C^) and YP"" = C.^^^^. 



(4.13) 



With this, the solution of equation (3.29) is absolutely straightforward. Since g^u\s = duui 
we have Sg^u = in the equation and if (4.6) and (4.7) hold, the following equalities must 
also hold on the boundary S: 



D^dg^s 



-2g<^5A% and 6{Dr,e)\s = C^^' 



TTj' 



(4.14) 



The equation for [7^°" reduces thus to a linear differential form in SA"^ . With (4.14) and 



^9fiu = 0, the equation for the superpotential (3.29) may be written 



SUf" 



d6^C 1 5£/f^r 



dB 



A 

fll/cr] 



3d{D,^gn^y 



rju 



6A 



fJbV 



(4.15) 



In this expression it is understood that the square brackets [p and a] mean that the expression 
is anti-symmetrized in pa. Since the factor of 5A^^ contains only the metric and the curvature 
tensor which on the boundary are given thanks to (4.5) and (4.6), provided the T's are 
derivable at least once, we may write that 



UP" 



d6^C 1 sr/f^r 



dB 



A 

fll/cr] 



3d{D,^gn^y 



rjU 



A 



fJbU' 



(4.16) 



s 



jpa 



-A 



Notice that U =0 because A = 0. U^" is the superpotential for theories of gravity 
derived from a Lovelock Lagrangian. 
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5 Application to General Relativity in N dimensional space- 
times 



(i) The variational derivatives 

Let 

B ^ g^'B^, = b^'^P^B^.p,, (5.1) 

where 

Then, if a hat Uke in g^^, means multiphcation by \/—g, the "relative" Einstein-Hilbert action 



isi4 



Ai = I Cid^x = ^ i{B- B)d^x. (5.3) 



The variation of the action, 



5 A, = JiC^^Sg^" + Cr^^\.)d''^ -If b^'\SA\^dS„ (5.4) 

The variational derivatives, 

^.^=^ = ^(%.)-^5.^^) and Cr = ^^=^y{D,g-^^-D,gn- (5-5) 

We may and have indeed used S^jy in terms of background derivatives and in terms of A's so 
as to keep every term covariant. The boundary conditions that would follow from the Action 



principle would be b'^^^^SA dSf^ = 0. Since neither b'^^'^ ^ = nor 6A dS^ = are, in 
general, compatible with (4.5) and (4.6), we shall add a divergence and take the following 



new action^^ 



A[ = A, + -Ji b^''\A\^dS,. (5.6) 

^ Js 



The variation of this action is 



SA[ = j\c,Jg^^ + CrdA^M""^ + \j 5{b^P\)A\jS^. (5.7) 

The boundary conditions may now be 5{u'^''^ = or b^'^P'^\s = g^^yPg''\^ and these conditions 
are compatible with (4.5). 



^''The 1/2k factor is for convenience. 

^^This is not the only surface term possible which leads to acceptable boundary conditions. We might also 



add \ /j, hf^P^xA^p^dS^ instead to Ai. 
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(a) The superpotential 

The superpotential (4.16) in this case is simple to calculate: 

Ur = -iH'''^\\,A\, = -pfe^^^^l^^^ff.A. (5.8) 



The last equality follows from (4.3). [/("^ is, as expected, the KBL superpotential. This may 
not be apparent because the KBL superpotential is more often written locally like this [20] : 

U^^ = - (d^p^^ - D^l'^l) + -i^fk"^ where k'' = g^^A"^^ - g'^'A'^^^. (5.9) 

To see that (5.9) is the same as (5.8) we notice that k'^ can be written as follows: 



j^a ^ 2b''^'\A^^^ and D^P^''^ - D^p^^ = {g^'^P - g''^P)D^C''^ + gt^^PA^l^^. (5.10) 

But since g^u\s = gui/i the first term on the right-hand side is zero. So, taking (5.9) into 
account and using the expression for the tensor b defined in (5.2), one easily finds that (5.10) 
is the same as (5.8). 

One may also write (5.8) like this: 



ur = k^D,{r^pr^''-r^prh. 



(5.11) 



It appears in this form, on a sphere at spatial infinity, in [1] and [9] with gP'^ to leading 
order in 1/r. The physical properties of the KBL superpotential in General Relativity are 
summarized in [18]. It is worth noting that the superpotential is valid on any S, whether at 
spatial or at null infinity. 

6 Application to Gauss-Bonnet theories in the Palatini rep- 
resentation 

(i) Lagrangian, variational derivatives and boundary conditions 

Gauss-Bonnet theories of gravity^^ in a Palatini formulation apply to spacetimes with 
more that 4 dimensions. Another difference with Einstein's theory of gravitation lies in the 



^^Gauss-Bonnet theories of gravity are rarely considered alone. What is usually used is Einstein's theory in 
N dimensions to which a Gauss-Bonnet term is added. This is then called an Einstein-Gauss-Bonnet theory 
which we consider briefly below. 
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relation between the metric and the connection coefficient. In Einstein's theory C^ = 
imphes D\g^^ = 0. In Gauss-Bonnet theories D\g^^ = is always a solution of L^" = 0. 

Equations (4.1) to (4.6) and (4.10) to (4.16) hold in Gauss-Bonnet theories. The varia- 
tional derivatives are derived from a Lovelock Lagrangian of order two which we took from 
Jacobson and Myers [16]. Define 

S= 9^''9'"'{B'^^xp^ „^u- B'^^^^B p^j^-2B p^^Byx + B^yBfj^- Bp^By„). (6.1) 

Then, the relative action with respect to the background, 

A2= I C2d^x= ^ l'{B-I)d^x. (6.2) 

The variation of the action is of this form: 

5A2 = JiC^uSg^'' + r/'JA^^Jd^x - i ^ P^'\SA^p,dSp; (6.3) 

in this the variational derivatives with respect to the metric components are as follows: 

^M^ = J^ = "5 ^^a^'^ + ^/^^^^u (6.4) 

in which 

2k£^^ = 9'""{B'^pX(,,^%va + B\:,)f,B^v + BpaB(^^^) + B'^^^^^B^)^ - Bp(^^B^y - B"^ ^^^^^B\^^) , 

(6.5) 
while variational derivatives with respect to the connection coefficients are of the form 
D[0{g)B] where O depends on the metric components only, the Bianchi identities elimi- 
nate the derivatives of i?'s and the result is of the form D[0{g)]B] explicitly: 

2k£/'^= 2k^ = lD^[g^'(P5lg'^^^ + g'^i-Slg^^PjB^.r + SDr,{gf'(-6lg^^(p)B%r 

'^ pa 

+ lDn{5[^b'''^^PP + 5^fhP'^>^')Bp^ (6.6) 

The boundary term introduces a tensor denoted by Davis [5] as P : 

ppupa^^ = {RP^P'' - 2g^^PR''^'' + 25'^t^i?'^]^ + Rg^^Pg''^''). (6.7) 

ppupa^ which has manifestly the symmetries of the curvature tensor like 6^'^'"^, is also diver- 
genceless: D^P'^'P'' = 0. 
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We may take a new action which is indeed compatible with (4.5) and (4.6) by adding a 
divergence to L^'- 

4 = / L2(fx + - / P^''\^^^„dS^,, (6.8) 



so that 



M'2 = jiC.xdg'^' + C^'''5/\\„)d''x + -^j> 6{P^P\)A\^dS^. (6.9) 

The boundary condition that follows from the variational principle is now 5P^^'^^" = which 
is equivalent to (6.7). This condition will automatically hold if (4.5) and (4.6) are satisfied 
and the T's are derivable at least once. 



(ii) The superpotential 

The calculation of the superpotential 1/2^^ given by equation (4.16) is somewhat more 
complicated than for Einstein's theory in N dimensions. MathTensor and Mathematica were 
valuable tools to check our calculations. A condensed but readable formula for t/^"^ is perhaps 
this one: 



{A^'^^uR'^Ki"' + R^'^'Ku e^ - ^''("A^, e")), (6.10) 



the factors of the A's are, of course, to be evaluated on the boundary, that is, in terms of the 
background geometry which is generally much simpler than that of the spacetime itself. 

Deser and Tekin [8], [9] have calculated the mass-energy in generic higher curvature 
gravity theories, in particular on anti-de Sitter backgrounds, motivated by the role of these 
backgrounds in string theory. 

Because of its importance in string theory and as an opportunity to relate our results 
with previous calculations, we consider now what becomes of C/j*^ if the curvature tensor of 
the background is of the form 



— 1 2 - 

Rfii/pcr = 'ToydnpQau ~ 9ii<j9pi/) ~ To^t^i^pcr- V"-11J 



j2\i>ppycn/ ypaiJpu) ~ j2"t^i^pcr 

With such a background (6.10) becomes rather simple:^^ 



U^^ = ^{N - 3){N - 4)C/f". (6.12) 



'^ Formulas (6.10) as well as (6.12) are at variance with the superpotential suggested in [7]. 
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We now consider the Einstein-Gauss-Bonnet theory of gravity and follow Deser and Tekin 
[9]. Our action is written like this 

A = hci + aC2)d^x + - I 2Ao(y=5 - y^)d^x. (6.13) 

a is a coupling constant and Aq contributes to the overall effective cosmological constant A 
which is related to"*^^ /^ as follows: 

,,_(iV-l)(iV-2) ^^^ 

2A 

The superpotential in this case is of the form 

In the case considered by Deser and Tekin, 

-2/2 
" = -"^^ = (iV-3)(iV-4)- ('•'') 

Thus, 

UP"" = -U^". (6.17) 

This is indeed their superpotential on a sphere at spatial infinity, as can be figured out from 
their formula (31) in which we must set a = /? = 0. 

7 A brief summary and some comments 

(i) First the summary 

The superpotential Un" of a Lovelock Lagrangian of order n, £„, n = 1, 2, 3, • • • , or the 
superpotential U^" = ^^ aJ^Un" of a linear combination of such Lagrangians C = '^^ a"/^™, 
with coupling constants a", is given by equation (4.16). The structure of Un'^ is rather 
obvious. Un'^ is homogeneous of order (n — 1) in the curvature tensor of the background and 
is linear homogeneous in A^^ = r^^. — T^^. 

We calculated Ul'^ which is well known and has been obtained by various authors with 
different methods. It was, however, useful to show the reader that we recovered at least well 
known results. One practical novelty, besides equation (4.16), is an explicit expression for Uq^ 



in [9], I is defined with the opposite sign so that their I for AdS spacetimes is negative. 
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on arbitrary backgrounds, formula (6.10). Another one is an expression for C/2"^, see (6.12), 
on anti-de Sitter backgrounds. This later expression provides the same superpotential as 
that found by Deser and Tekin at spatial infinity, which they obtained using a very different 
method. The result gave us further confidence that the procedure of Silva, extended in Section 
3 to Palatini's representation of gravity fields, did indeed work and that equation (4.16) was 
correct. One might have calculated f/g"^ and U^'^ using £3 and £4 explicitly written in an 
Appendix of a paper by J. T. Wheeler [32] ^^. This would have been hard work with little direct 
prospect of applicability. The method of Section 3 may of course be applied to variational 
derivatives of higher order that are not Lovelock Lagrangians, the type of Lagrangians that 
interest string theorists. 

(a) Now some comments 

a] Regarding [/("^ as given by equation (5.8). On a flat background, in Minkowski coordi- 
nates and with ^ the Killing vector of spacetime translations, [/^'^ is exactly the superpo- 
tential that Freud [31] found, almost 70 years ago, to calculate mass-energy and total linear 
momentum. One wonders why Freud did not calculate the angular momentum on the same 
occasion. 

b] We emphasized several times that the superpotential holds at null as well as at spatial 
infinity and that [/{"^ at null infinity gives the Bondi mass. We do not know if anybody 
got interested in radiating fields for Einstein-Gauss-Bonnet theories. C/2"^ is, of course, the 
superpotential appropriate for calculating the Bondi mass. 

c] We imposed Dirichlet boundary conditions. Julia and Silva showed that imposing Neumann 
boundary conditions lead to Komar's superpotential^*^. Neumann boundary conditions have 
been considered in recent works. See for instance a paper by Kofinas and Olea [23] on 
Lovelock anti-de Sitter gravity. Our superpotential 1/2'^ does not apply to such spacetimes. 

d] Jacobson and Myers used the first law of thermodynamics to define mass energy of Lovelock 
black holes. Mass-energy has, however, little to do with thermodynamics because asymptotic 
spacetimes ignore the source of gravity. A direct calculation of the mass-energy of Lovelock 



^^See also Miiller-Hoissen [25] who gives £3. 



^"A generalization of Komar's superpotential, that does not need a background, to any type of boundary 
conditions, can be found in several papers of Obukhov and Rubilar's ; see for instance [26]. The role of 
backgrounds is here replaced by another ingredient, "generalized" Lie derivatives. 
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black-holes, independent of thermodynamic considerations, is naturally provided by superpo- 
tential JJP" . An approach, similar to that of Jacobson and Myers, was used by Gibbons et al 
[14] to find the mass of Kerr-anti-de Sitter black holes. Needless to say the same mass-energy 
is obtained with the KBL superpotential C/f^, see [6]. 

d] Finally, we are still wondering about the following problem. Noether's identity provides 
a conserved current from a given Lagrangian. The conserved current is the divergence of a 
non-unique superpotential. The form of that current (see Section 1) tells us what physical 
meaning to attribute to closed surface integrals of the superpotential assuming it has been 
properly chosen. With Silva's prescription, things work the other way round. One calculates 
a unique superpotential but one must look at its divergence to figure out the physical meaning 
of its total flux. Both methods are sound and well defined but what is the formal connexion? 
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